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Abstract 

Quantum mechanics sets severe limits on the sensitivity and required circulat- 
ing energy in traditional free-mass gravitational- wave antennas. One possible 
way to avoid these restrictions is the use of intracavity QND measurements. 
We analyze a new QND observable, which possesses a number of features that 
make it a promising candidate for such measurements and propose a practical 
scheme for the realization of this measurement. In combination with an ad- 
vanced coordinate meter, this scheme makes it possible to lower substantially 
the requirements on the circulating power. 

I. INTRODUCTION 

In we presented an analysis of two qualitatively new schemes for the extraction 

of information from free-mass gravitational- wave antennas Common features of these 
schemes are the use of nonclassical quantum states of the optical field inside the resonators 
and of QND methods for intracavity measurements of the variations of these states. This 
becomes possible only with the realization of optical field relaxation times r* much longer 
than the measurement time Tmeas — 10~^ -7- 10~'^s. One significant advantage of intracavity 
measurements is that they require lower levels of circulating power than traditional schemes 



with an antenna with a coherent pump. In |T[ and in the subsequent article by Levin [Q, 



easily, was to place an additional partially transparent mirror-probe mass at the intersection 
of the two arms of a gravitational antenna. This results in the formation of two coupled 
Fabry-Perot resonators. Displacement of the end mirrors under the action of a gravitational 
wave leads to a redistribution of the energies in the arms, which pushes the central mass. The 
absolute displacement under optimal conditions is simply equal to the relative displacements 
of the end mirrors {hL/2, where L is the arm length and h is the amplitude of the variation of 
the metric), and the light in the system behaves like a rigid bar. The displacement associated 
with an independent mass that does not interact with the optical field can be registered 
without consuming a large amount of power. A rigorous general relativistic justification of 
the schemes in can be found in ||5[. 

The merits of this intracavity measurement are the following: a) in the resonators, the 
required nonclassical quantum state (close to a Fock state) is formed automatically; b) direct 
measurement of a displacement hL/2 consumes relatively little power; c) precision higher 
than the standard quantum limit can be obtained. 

In 0, we did not make an analysis of the minimal energy of the optical field S in 
the system required to preserve the sensitivity. Another important unanalysed problem is 
the connection between the achievable resolution and a chosen procedure for displacement 
measurement. 

It is important to note here that the provision of substantial values of £^ is a key problem 
for large-scale gravitational wave antennas, and that this problem has not a technical but a 
fundamental nature. Indeed, the proposed sensitivity levels of such antennas will be close 
to the standard quantum limit for the displacement of the masses M of the end mirrors: 



where Ugr is the frequency of the gravitational signal and Tgr is its duration (we omit in 
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the optical cubic nonlinearity x^^^ of thin plates inserted in an antenna was exploited. 

The idea of our second scheme , which, in our opinion, can be implemented relatively 




n 



(1) 



our estimates numerical terms of the order of unity that depend on the form of the signal). 
According to the Heisenberg uncertainty relation, the momentum should be perturbed by a 
value of the order of: 

h 

2AXSQL 



= J^Mu^^Tgr (2) 



This perturbation must be provided by the uncertainty in the energy S in the interferometer, 
which, thus, cannot be less than 



AS = .r^LAp = L^hMul, (3) 

This value is not especially large; for example, for L = 4 x lO^cm, ujgr = 10^s~^, and 
M = lO'^g (the parameters of the LIGO antenna), 

A£ ^ 4 X lO'^erg, (4) 

and in the case of nonclassical states of the optical field, in which AS ~ S, the necessary 
resolution can be obtained at very low energies. However, for coherent states in which 



AS = ^JhuoS, (5) 
where cUo is the optical frequency, the requirements are very strict: 

£ . ^i^. (6) 
For the same parameters as before and Uo = 2 x 10^^s~^, 

SsQL ^ W'erg, (7) 

and if Ugr = 10^s~\ then 

SsQL ^ lO^'erg. (8) 

In this paper, we analyze a new intracavity scheme that is, in some sense, complementary 
to the "optical bars" scheme. In this scheme, the optical field forms in a quantum state that 
is close to states with squeezed phase; this is known to allow, in principle, a dramatic 
decrease in the optical quanta because S^p ~ (Non-QND measurement of a similar 

observable was proposed in |P]). 



II. A CROSSQUADRATURE QUANTUM OBSERVABLE AND A SCHEME FOR 

ITS MEASUREMENT 



The basic idea of the new scheme for an intracavity readout system is the use of two 
modes excited in the Fabry-Perot resonators of the antenna's orthogonal arms. If the modes 
are not hnearly coupled (this is critical in this scheme), they can be tuned as close to each 
other as {ui —uj2)Tmeas ^1- As a result, the frequency variation in one (or both) resonators 
produced by a gravitational wave will lead to the appearance of a phase difference with the 
oscillation amplitude 

5^ ^ ^, (9) 

which we propose to register. Since no meter has been invented thus far to directly regis- 
ter the phase difference between two quantum electromagnetic oscillators, another variable 
proportional to 5ip is required. 

We propose to measure the averaged product of the two quadrature components of two 
different oscillators, which, in the limit of large numbers of quanta, is very close to a phase 
measurement. One possible scheme for the realization of the proposed crossquadrature ob- 
servable is depicted in Fig.l. This scheme is based on the use of ponderomotive nonlinearity 
in a way similar to that in |^. Mirrors A' and B' direct the optical beams reflected from 
the end mirrors A and B and transmitted by the 50% beamsplitter C on opposite sides 
of the double highly reflecting (zero transmission) mirror D (to eliminate linear coupling). 
In the engineering realization of this scheme. A' and B' can be rigidly connected to the 
beamsplitter, and can be focusing reflectors, making it possible for the mass m of D to be 
smaller. 

It is easy to see that, due to the beamsplitter, the optical beams from arms A-C and B-C 
interfere in the shorter arms such that one of them has amplitude proportional to ai + ia2 
and the other has amplitude proportional to 02 + iai (ai^2 are the complex field amplitudes 
in the longer arms). This is valid if the geometrical conditions in Fig.l are satisfied. As a 
result, the ponderomotive force Fpond acting on mirror D will be proportional to: 



Fpond oc \ai + ia2f — \a2 + iai\'^ ~ 4|ai| |a2|5(/7. (10) 

Provided that the initial optical energy S/2m the two arms is nearly the same {S — fiuJoN — 
hcUoCi^, a = |ai| = |a2|), in a quasistatic approximation, this force will be 

g 

Fpond ^ jScp (11) 

Note here that there is no direct linear coupling between modes in this scheme. In other 
words, modes in the resonator are coupled via the x^^^ nonlinearity resulting from the pon- 
deromotive effect. Linear coupling is due only to the movement of the mirror D. The shift 
of D changes the lengths of the shorter arms, changing the interference conditions on the 
beamsplitter, which consequently leads to a redistribution of the optical photons between 
the two modes. 

This scheme reahzes indirect QND measurement of the operator 

Xn/2^i{ata2- atai), (12) 

where 0,12 and ai 2 are the creation and annihilation operators for two different oscillators 
with the same frequencies uu. The operator X^/2 presents a special case of the family of 
operators 

Xe = 0^026^^ + atciie-'^, (13) 

which we propose to name crossquadrature operators. These operators commute with the 
Hamiltonian of the two modes: 

[Xe, huj{a^ai + 0^02)] = 0, (14) 

i.e., they are, indeed, QND variables. The eigenstates of the crossquadrature operators have 
the form 

\N,n) = , ^ (a+ + a+e-^Y(Qi + a2e'^)'^"lO). (15) 

^j2^n\{N -n)\ 
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where N is the sum of quanta in the system and n is an integer in the range from to 
A^. In this state, each of the quanta has equal probabihty to reside in either arm of the 
interferometer. However, the amphtudes of these probabihties for n quanta are orthogonal 
to those of the other N — n quanta. Due to this peculiar entanglement between the modes, 
we shall call eigenstates of the crossquadrature operator symphotonic quantum states. 

The eigenvalues of the operator Xg are n — (A^ — n) = 2n — N, i.e., measuring the 
crossquadrature variable, the observer determines the difference between the two kinds of 
quanta. Symphotonic states ([T5|) are very sensitive to the change of the phase difference in 
the two oscillators. As we show in Appendix A, a phase shift leads to a transition between 
states with different n (preserving the total number of quanta), that can be detected by 
measuring The probability of this transition is equal in the case 6ip <^ 1 to 

p = ^{N + 2n{N -n)), (16) 

and when n ~ N/2, p tends to unity when S^p ^ y/S/N, thus allowing, in principle, the 
registration of these small phase shifts. 



III. LIMITATIONS ON THE SENSITIVITY 

It is not difficult to show that the finite masses of the mirrors A' and B', as well as the 
mass of the beamsplitter C, do not influence the behavior of the system if these masses 
are substantially greater than the mass m. We will use a standard linear approximation, 
in which the optical field can be represented as the sum of the large classical dimensionless 
amplitude A and the quantum annihilation operator a, neglecting terms of the order of 
and higher. We suppose also that r* and relaxation time of the mass m is very large in 
comparision with other characteristic times. In this case, the equations of motion will have 
the form: 

ddi{t) fixi{t)—x{t) ih{t)\ f i5o 



dt 
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da2{t) 
dt 



^ ~ir (""'^ ^ + ""^^ ^) 

^—d^ = + at{t)) (17) 

where Xi_2 are the displacements of the mirrors A and x is the displacement of I?, 5o — 
1/2t* is the decrement of the optical losses in the resonators, 6i,2(<^) are the corresponding 
annihilation operators for the heatbath modes, which satisfy the commutational relations 

6^2(^0] =<5(^-^'), (18) 

pmeter jg ^j^g fluctuational reaction of the coordinate meter on the mirror D with mass 
m, h{t)/2 is the relative change of the optical lengths of the resonators (in the case of a 
gravitational antenna, this is the dimensionless metric variation), and F'^^'^^ is the Nyquist 
fluctuational force acting on the mass m. 

The characteristic equation of this system is: 

+ = 0, (19) 

where 

-f^V'°. (20) 

It has roots with positive real parts of the order of v. Thus, there exists in the system a 
dynamic instability with a characteristic time To suppress this with a feedback loop, 
it is necessary to have 

U < UJgr- (21) 
The signal-to- noise ratio is equal to (see Appendix B) : 



n J m^{u^-uj^)^S^ + 27rvjj'^{iy^-uj^)Sa:F + uj^{SF + S,n + So) 27r'' 



where huiiuj) is the signal spectrum, 



Sm = - (23) 



is the spectral density of F"^ {k is the Boltzmann constant and T is the temperature), 

is the spectral density of the fiuctuational force due to dissipation in the optical resonators, 
Sx and Sp are the spectral densities of the additive noise x'"^*'^^ of the meter and of i^'"^*'^^^ 
and SxF is the cross spectral density of x"^^^*^"^ and F^<^^^^, The values of Sp, and Sxf{^^) 
must obey the Heisenberg inequality 0]: 

Sx{u;)Sp{u;)-SU^)>^. (25) 

The condition for the detection of a signal can be represented in the form: 



h > yh'^neter + ^mec/i + ^opt: (26) 

where 



_ LUgr 2KTm _ Ugr / 2kT , . 

iT-mech — \ — ~ '^~a'\ ~ ^7 l^'J 

^^o V ^mV ^ V ^^V^ 



is the limit due to the thermal noise of the mass m, 



is the limit due to the optical losses (r^j. is the duration of the signal), and hmeter is the limit 
due to the quantum noise of the meter. 

It is important to note that the limitation ( pHj) is also valid for the previous scheme 0, 
based on a different principle for intracavity measurement (this follows from formula (10) of 

§)• 

The value of hmeter is determined by the magnitudes of the spectral densities 
Sxiyj), Sf{uj), and Sxf{^) and their frequency dependence. In the case of a plain coordinate 
meter: 
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Sx{u!) = const, Sf{uj) = const, Sxf{^^) = 0. (29) 
With regard to limitation (pT]), values corresponding to optimal tuning of the meter will 



be: 

Sf = Sx = -. (30 

The ultimate sensitivity of the meter is determined, in this case, by the formula 

L 



h 



meter 

to, 



^{^) hsQL{M). (31) 



Thus, because of it is impossible in this case to reach a sensitivity corresponding 

to hsQL{M). 

To preserve a sensitivity at the level of hsgiiM) and lower the requirements on the 
energy, one can use an advanced meter providing higher precision for monitoring the mass 
m. A speed meter can be used for this monitoring. This can be realized in the form of 
an ordinary parametric electromagnetic displacement transducer (operating at microwave 
wavelengths) with an additional buffering cavity, coupled with the main (working) cavity 
We show in Appendix C that, in this case, 

•^^•^ = A 2 u/ ■ Sf{uj) = , Sxf{uj) = -- cot $, (32) 

where Ug is the microwave frequency. We is the microwave pump power, d is an equivalent 
parameter with the dimensions of length, which characterizes the tunability of the transducer 

^ =^9^' ^^^^ 

Qe is the beat frequency between the working and buffering resonators, which must satisfy 
the conditions Qe ^ ^gr and ^e/T* ^ ^gr (K the relaxation time due to the coupling 
with the transmission line), and $ is the phase of the local oscillator used for detection of 
the microwave signal. For optimal tuning of the meter parameters, when 
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and 



= a^^L (34) 



cot$ = -^, (35) 



the limiting sensitivity will be 



h^eter = V2hsQL{M). (36) 



If, for example, = 2 x lO^^s^^, L = 4 x lO^cm, cUg^ = lO^s^^ (these values correspond 
to the values for the LIGO antenna 0), m = Ig, and £^ = lO^erg, then ~ 5 x 10^s~^, and 
condition (^) is satisfied. If in addition d = 1cm (the value achieved in high-Q sapphire 
disk resonators |jlll) and fig = 3 x lO'^s"^, then the required microwave pump power will be 
We = 3 X lO'^erg/s. Thus, the analyzed scheme makes it possible to dramatically decrease 
the requirements for the optical circulating energy by using a microwave transducer with a 
reasonable set of parameters. 

Under these conditions, however, the requirements for the level of dissipation in the 
probe mass m increase as the signal that must be registered decreases: 

'^Vm ^ —AVSQL, (37) 

where 



^VSQL 



\ 



* (38) 



rriTgr 



If the above parameters are chosen, Avm — 1/8Avsql- In order for the dissipation not to 
deteriorate the sensitivity, it is necessary that hmech < hmeter, or 

< ^ir- (39) 

m 

For example, for T = AK, > 3 x lO^s. Thus, the requirements for the dissipation in the 



mass m are severe, but achievable 12 



10 



Losses in the optical resonator will not influence the sensitivity if hopt < hmeter, which is 
equivalent to the condition 



or, for the parameter values introduced above, r* > Is. This is quite possible in a LIGO-type 
detector with optical mirrors available today. 



In 0, not all regimes for the "optical bar" scheme were analyzed in detail. Moreover, 
there was unfortunately an error in the formula following formula (12) (term "1" under the 
root should be omitted). 

Here, we shall limit our treatment to "wideband" regimes, when the range of the signal 
frequencies is far from the resonant frequencies in the signal-to-noise integral; these regimes 
are the most useful from the practical point of view. The "narrowband" regime, in which it 
is possible to attain sensitivity better than the SQL, has already been considered in detail 
in [0. In our analysis, we shall assume that ujgr < Q {Q is the beat frequency in the system 
of two coupled optical resonators; the case of > f2 is difficult to realize in practice, and 
does not provide any interesting new results) and m <^ M. The behavior of the system is 
determined by the parameter with the dimensions of frequency 



This frequency describes the influence of the ponderomotive force on the dynamics of the 



It is possible to distinguish three cases, depending on the level of the circulating energy 
(the value of 9). 



IV. COMPARISON WITH THE "OPTICAL BAR" SCHEME 




(41) 



system, and plays a role analogous to u (see formula (PP])). 



A. Weak pump power, < UgrO, 



If a plain coordinate meter is used, the calculations give the following result: 
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IOmeter 



(^) ^SQLim) > hsQLim), (42) 



where 



hsQL{m) = ^ 



n 



(43) 



If a speed meter is used, the sensitivity can be higher: 



''meter 



02 

but is still lower than even hsQhijn)- 



hsQL{m) > hsQLim), (44) 



B. Intermediate case, ujgrQ < < ujgr^^2M /m 



For a plain coordinate sensor, the best sensitivity in this case is 



hmeter = -^hsQlim), (45) 



UJgj.^ . 

i.e., hmeter IS Smaller than hsQiija), and hmeter hsgiiM) if 6^ ujgrQ^2M/m. The 
required optical energy in this case is 

e ^ ^Ss,. > S.,.. (46) 

Ulgr 

The use of a speed meter in this regime does not give a gain in sensitivity, however an 
increase in sensitivity is possible if an advanced coordinate detector with correlated noises 
is used {SxF 7^ 0) In this case 

hmeter = (^^^^^) hsQlim), (47) 



if 92 < ujgrVt^2M/m, and 

hmeter = hsQL{M) (48) 

otherwise. The required energy in the latter case can be lower than Ssql, but with respect 
to a possible dynamical instability, which appears when 0^ > f2^/4: 
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C. Strong pump power, 6^ S> ujgr^^2M/: 



m 



This is the "optical bar" regime, when the masses M and m move together, and are 
connected by electromagnetic rigidity. In this plain coordinate meter provides a 

sensitivity corresponding to the standard quantum limit hsQL{M), and use of a speed meter 
makes it possible to overcome this limit, but with higher energy: 



Note that, in this case, also, the total mass ~ 2M is present in the expression for the 
thermal limit. This result is quite understandable, since, in this regime, thermal fluctuations 
of the small mass m act on the large compound mass 2M + m. 



Quantum mechanics sets severe limits on the sensitivity and the required circulating 
energy in traditional free-mass gravitational-wave antennas. One possible way to beat these 
limits is to use intracavity QND measurements. In this paper, we have analyzed a new QND 
observable and its corresponding symphotonic quantum states, which possess a number 
of features that make it promising for experiments requiring registration of small phase 
variations: 

1) Unhke other known QND observables, this one is a joint integral of motion for two 
quantum oscillators with equal frequencies. 

2) The crossquadrature observable is very sensitive to the phase difference of the oscil- 
lators. Phase differences of the order oll/N (the theoretical limit for phase measurements) 
can be detected, where N is number of quanta in the system. 

3) Well-known methods for the QND measurement of electromagnetic energy can be 
used to measure this new observable. 

We have considered a practical optical scheme in which the new observable can be used 
for the detection of gravitational waves. Our estimates show that, in combination with 




(50) 



V. CONCLUSION 
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advanced coordinate meters, this scheme provides a sensitivity of the same order as that for 
planned antennas at significantly lower energies. 

Summarizing the results of this article and of we conclude that intracavity mea- 

surements with automatically organizing nonclassical optical quantum states make it possi- 
ble, in principle, to lower the required power levels and in several cases to achieve sensitivity 
better than the standard quantum limit. 

We note also that the schemes we have analyzed do not cover all possible geometries for 
intracavity measurements with ponderomotive nonlinearity. Better realizations with higher 
responses are probably possible. 
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APPENDIX A: THE EVOLUTION OF A SYMPHOTONIC STATE 

The evolution operator describing the phase shifts in the first mode and 02 in the 
second one is equal to 

^^(01,02) = exp I — 1 (Al) 

where rii 2 are the operators for the number of quanta in the modes. Hence 

U{h, h)atU^{<l)u 02) = ate-"^\ (A2) 

?7(0i, 02)a+f>+(0i, h) = a+e-*^^ (A3) 



and 



f/(0i,02)|O) = |O). (A4) 



Taking into account formula ([Iq) and omitting the unimportant factor e 2 , we can 
obtain: 

t/(0i,02)|iV,n) = ^ =(i+cos(50 + zE+sin(50)"(^+cos(50 + iA+sin(50)^-"|O) 

(A5) 

where 

i+ = a+ + a+e-^^ = a+ - a+e^^^ (A6) 

and 



■>2 - VI- 



(AT) 



If (50 < 1 then 



f/(0i,02)|Ar,n) ^ (^-^-^{N + 2n{N -n))^ \N,n) + 
i5(t) I^^Jn{N -n + l)\N, n - 1) + ^ {n + 1){N ~ n)\N , n + l) 
^ - 1)(A^ -n + 1){N - n + 2)|A^, n-2) + 



^{n + l){n + 2){N - n){N - n - l)\N,n + 2)j (A8) 
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Thus, the probabihty for changing the number n is equal to 

p = 1 - |(iV, n\U{<Pi, 02)|A^,n)P ~ ^(A^ + 2n{N - n)). (A9) 



APPENDIX B: SIGNAL-TO-NOISE RATIO 

Equations ([T7|) can be rewritten in the form: 

){t)+A I {{btH - &^(cu))e'("+"°)* + h.c.) du 



— -X ' ^ ' ■ ' 



dt L 

^ y{t) + uJoNh{t) +aJ tJ'^ {{btiu) - 6+(cu))e*("+"°)* + h.c.) duj 











d'^x{t) flUJo ^ fpmeter(.\ 

rd'^yjt) _ huo 

dt^ 



= (Bl) 



where 

Af = A{di + df — d2 — a^), = iA{d2 — — cii + di), (B2) 

A/" - difference of number of quanta in the two arms, y = [xi — X2)/2, and h.c. stands for 
Hermitian conjugation. Hence, the spectrum of x{t) is equal to 

X{uj) = Xsignali^^) + X (u) , (B3) 



where 



hujiN —ILO' 



3 



is the signal spectrum, 

^^'^'^ ^ m{v^ - uj^) (^"'*"(^) + + F°P\oo)) (B5) 

is the spectrum of fluctuations of x{t), and F°'p^{ijj) is the spectrum of force 







(B6) 
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The output signal of the coordinate meter is equal to 

= x(t) (B7) 
where a;'"^*^^(t) is the additive noise of the meter. Hence, 

X{UJ) = Xico) + X^^'^^UJ) = XsignaM + X{uj) + x'^^'^^u), (B8) 



and the SNR is equal to (|2|). 



APPENDIX C: MICROWAVE SPEED METER 

Let us consider two coupled microwave resonators. The first is connected to an output 
waveguide (see Fig.2), and it's eigenfrequency depends on the coordinate x to be measured: 

^{x) = c., (l - ^) , (CI) 

while the second is pumped by the resonant power Uq cosuet. 
The equations of motion for such a system are 

^^^^ + Ujlq2{t) = 2cUefiegi(t) + — Uq COS UJet, (C2) 

ac^ p 
where gi 2 are the generalized coordinates of the resonators, p is the wave impedance of the 

resonators, 6e = 1/2t*, t* is the relaxation time of loaded first resonator, and U^^"^^ are 

fluctuations in the waveguide (we neglect intrinsic losses and corresponding fluctuations of 

the resonators). 

Linearizing these equations in the strong-pumping approximation and using the method 
of slowly varying amplitudes, we can obtain: 

— h deai[t) = -ilj}2[t) 

at p 

dhiit) , , , , ujeqo , , /X fZ/'^'^Ht) 



dt d 
da2(t) 



dt 
dhit) 



dt 
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-^ebi{t) 

n^ai{t) (C3) 



where ai,2 and 61,2 are the amphtudes of the cosine and sine quadrature components of gi,2, 
f//^""^* are the same for f/-/^'"^*, and go is the mean value of the amphtude of oscillations in 
the first resonator. 

Solution of these equations in the spectral representation gives: 

where C^u) = Ql — + iuSe- The output wave in the waveguide can be represented in the 
form: 

jjout^^-^ _ jjfiuct^^-^ _ 2(5eP dqi(t) _ 

UJe dt 

(?7f ""*(t) - 25epbi{t)) cosiOet + + 26epai{t)) sincu^t (C5) 

If a homodine detector with Ulq sin{uet + is used, where $ is the phase of the local 
oscillator, the output signal of the detector is proportional to 

U{t) = (f/f - 25ePh{t)) sin $ + (f/f + 25epai(t)) cos $ (C6) 



Substitution into this expression of the solution ( p^ gives that the spectrum of U is equal 
to 

L[uj)a ^ ' 

where 

^meter^^^ = __ J __ (^^(^2 _ ^2 _ siu $ + [/f "'^*(cu) COS $ 



2ici;a;e(5egoPsin $ 



is the spectrum of the additive noise of the meter. 

The fiuctuational reaction force of the meter is equal to 



d 

or, in spectral form. 



pmeter^^^ = ^2^ai(t), (C9) 



^ ^'^^ IWd ^ ^ 

If the frequency u is relatively small: 

uj^^nl and uj6e<^nl (Cll) 

then expressions ( |C1U| , |C5D directly give the spectral densities 
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FIGURES 



A 



L-l 




FIG. 1. The scheme of measurement of the crossquadrature observable 




FIG. 2. The scheme of microwave speedmeter 
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Figure captions 



Figure 1. The scheme of measurement of the crossquadrature observable 



Figure 2. The scheme of microwave speedmeter 
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